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Abstract. The characteristic anti-linear (parity/time reversal, PT) symme- 
try of non-Hermitian Hamiltonians with real energies is presented as a source of 
i two new forms of solvability of Schrodinger's bound-state problems. In detail 

' we describe (1) their very specific semi-exact solvability (SES) and (2) their 

innovated variational tractability. SES technicalities are discussed via charged 
oscillator example. In a broader context, speculations are added concerning 
^ ■ possible relationship between PT symmetry, solvability and superintegrability. 
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■ PT symmetric quantum mechanics: A brief introduction 

Prehistory: Isospectral operators. Origins of the popular Bender's and 
Boettcher's PT symmetric quantum mechanics BO] IB 2) lie in perturbation theory 
[CGI IFGj . For an elementary illustration one may recollect the pioneering pa- 
per by Buslaev and Grecchi BG who proved the isospectrality of the Hermitian, 
_C ' spherically symmetric D— dimensional perturbed harmonic oscillator 

Is 

, (in its to— th partial-wave projection at any to = 0, 1, . . .) with its non-Hermitian 

' "unstable" anharmonic partner(s) in one dimension, 




H (UAO) = _d + ^ _ 2 _ + ^ _ _ > Q 

This means that the physics (and, in particular, the reality of energies) remains 
unchanged while the mathematics itself is significantly simplified when one weak- 
ens the Hermiticity of H^ PH °^ (g) = \H ( PHO ) (g)] ^ to the mere PT symmetry of 
its spectrally equivalent partners jDT) . The latter operators commute with the 
product PT of parity and time reversal. In the modern language of review M , 
one should rather speak about the pseudo- Hermiticity defined by the relation 



# AO »(j)=P W UAO \g) 



t 



P. 
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Recent progress: Exactly solvable examples. Paper |BG| extends the 
proof to the isospectrality between H^ PHO \g) and H^ UAO \g) to the imaginary 
couplings g = ih. The spectrum of the resulting self-adjoint double- well H^ UA °^ (ih) 
in one dimension (famous for its perceivably hindered perturbative tractability) 
may be then deduced from the "unstable" modification H^ PHO ' (ih) of the central 
quartic oscillator in D dimensions. The latter, non-Hermitian PT symmetric model 
is sufficiently simple in its ordinary differential "radial" representation of ref. BG , 

ar z r z 

r = r(x)=x — ie, s e (-00,00), h > 0, e > 0. 
In the context of the present proceedings this operator admits an immediate re- 
interpretation as a specific perturbed and regularized form of the most common 
superintegrable oscillator Hq UR °^ (0) of Smorodinnsky and Winternitz (SW, |FML 
[E|) in any number of dimensions and after its separation in cartesian coordinates. 

The latter observation is inspiring since it is known that the concepts of (max- 
imal) superintegrability and (exact) solvability are closely related at e — |T0| . 
In such a setting the weakening e > of the Hermiticity extends the class of the 
solvable models ZOj IZ1| as well as their interpretation in terms of Lie algebras 
|LC| and evokes a number of new open questions. We have addressed the first few 
of them in ref. |Z2j where we derived the exact spectrum and wave functions for 
the complexified SW oscillator H^ URO \o). Our next move in this direction was 
devoted to the closely related complexified (though still separable) Calogero model 
of three particles jZMj where an overlap of the integrability with PT symmetry 
acted as a generator of the new solvable self-adjoint Hamiltonians ZT . This expe- 
rience makes the further study of PT symmetric models of (super)integrable type 
very promising. In the present paper, we are going to describe some new results 
achieved in this direction. 

Running project: Partially solvable models. In place of the (not so easily 
feasible) analysis of the traditional quartic perturbations, we shall pay attention to 
the slightly simpler charged oscillators characterized, first of all, by their incom- 
pletely or quasi-exactly solvable 1 (QES) status. From the historical perspective, 
the acceptance of this concept was comparatively dramatic. In a prelude, vari- 
ous potentials have been shown partially solvable in terms of elementary functions. 
Thirty years ago, this period has been initiated by a two-page remark by Andre 
Hautot |H] who noticed that the charged oscillator possesses arbitrary finite mul- 
tiplets of exact Sturmian solutions in two (D = 2) and three (D = 3) dimensions. 
The Hautot 's choice of the elementary wave- function ansatz 

N 

<p(r)=e- r2/2 - 9r h nr n+K 
n=0 

(with arbitrary TV) preceded the discoveries of the QES sextic oscillator |SB| . of 
non-polynomial QES anharmonicities WF etc. During the "golden age" of the 
development of the subject, the existence of the common Lie-algebraic background 
of all the QES systems has been emphasized |T1| . A summary of the "state of the 
art" up to the early nineties has been offered by Alex Ushveridzc in his monograph 



^term coined by A. Turbiner, meaning an elementary solvability for finite multiplets of states 
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[Uj where the quartic polynomial oscillator is mentioned as a typical system without 
quasi-exact solvability. 

A change of the approach to QES models has been initiated by Bender and 
Boettcher |B1| who revealed and described the QES solutions for the "unstable" 
quartic polynomial oscillators. Their construction has been extended to all the 
partial waves in ref. |Z3j . In paper |Z4| we made the next move and turned 
attention to PT symmetrization of the Hautot's charged oscillator (cf. Appendix 
A). In the present continuation of this effort, we shall complete the picture by its 
extension, i.a., to all the partial waves (cf. Appendix B). In section 1 we describe 
our main result, viz., the quasi-even solutions which remained unnoticed in |Z4| . 
Unexpectedly, these "facilitated QES" states are available in infinite multiplets , all 
the elements of which are, basically, non-numerical. In a way emphasized in section 
2, these sets may easily serve as certain non-standard bases in some "sufficiently 
large" subspaces of Hilbert space. For this reason, we suggest to call them semi- 
exactly solvable (SES, cf. Table 1). 



Table 1. Tentative classification of solvability. 



class 


quasi — exact 


semi — exact 


exact 


solutions available 
range of couplings 
illustrative example 


finite set 
restricted 
H^ BBO \a,b) 


infinitely many 
restricted 
see below 


all 

any 
H (URO) (0) 



1. Quasi-exact solvability on complex contours 

1.1. Three-term recurrences. Let us consider the Schrodinger equation 



(1.1) 



dx 2 



\, ' + I —, — r + 2 I br{x) + r 2 {x) 
r^(x) r(x) 



tp n (x) = E n ip n (x) 



£ = £(L) = (L-l)/2, 



L = D - 2, D, D + 2, 



It contrast to the Hautot's QES problem of Appendix A it works with the purely 
imaginary charge and with a constant complex shift of coordinates r = r(x) — x—ie 
(cf. pQ). It also generalizes the £ = problem of ref. Z4 so that we must employ 
the more powerful ansatz 



N 



(1.2) 



^n(r) = e 



_ -r 2 /2-ibr 



n=0 



'such a feature is much more characteristic for the completely solvable models 
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After its insertion, the differential form of our PT symmetric bound state problem 
is replaced by the finite-dimensional matrix equation, 



(1.3) 



/ -S Co 
Ai B x 



V 



\ 



An-i Bn-i Cn-i 



A 



N 



with the real matrix elements, 

A n = b 2 + 2n - L - E, 
C n = (n + l)(n+l-L), 



B N 
A N+1 J 



B„ — 
L = 2l 



( Po\ 

P2 
\PN j 



= 



N < oo 



-(2n- 
I- 1, 



1-L)b- F, 

71 = 0,1,.- 



The last line of l|1.3f) represents a separate condition A^ +1 = which gives our 
energies in closed form, 



(1.4) 



E = E 



N 



2N + 2 — L + b 



N = 0,1,. 



The key consequence of this easy but important simplification lies in the emergence 
of a zero Cl-i = in the upper diagonal of the square matrix of the simplified 
system <|1.3fl ■ In Appendix B the presence of this zero enables us to construct the 
quasi-odd solutions at all I. 



1.2. The new, quasi-even solutions. All the coefficients in the polynomial 
wave functions l|1.2|l may be defined by closed formula (|4.2|) displayed in Appen- 
dix B. Let us now turn our attention to the quasi-even states. The superscript ( - + - ) 
will be introduced to mark their even quasi-parity which may be characterized, for 
our present purposes, by the simple negation of the odd quasi-parity condition l|4.3|l . 



(1.5) 



bo +) | 



\A +) \ 



\p ( l\\>0- 



For quasi-even states, the vanishing matrix element = plays a different 

role. It separates again the two subsets of equations but the upper one ceases to 
be trivial. This means that the related subdeterminant must vanish, 



(1.6) 



S M = det 



/ Bo Co 
Ai B x 



V 



A L -2 



\ 



Bl-2 Cl-2 

A L -i B L -i J 



= 0. 



Up to the exceptional, degenerate cases where eqs. (11.611 and (j4.4|l hold at the same 
time, we may drop the non- vanishing factor S^' 1 =/= and determine all the charges 



F 



A+) 



N k as roots of polynomial . 
The ^'—superscripted wave function coefficients are determined by eq. I|4. 2\ . 
this time in the full range of indices j = N, N — 1, . . . , 1. As long as the dimension 
of these determinants grows with the difference N —j, we shall recommend a return 
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to the recurrences 

/ -2N Pi - F 4-2L 
2 - 27V 13 2 -F 



V 



9-3L 

-4 /3iv-i 



\ ( P ° +) ^ 
p[ +) 



N 2 - NL 
Pn-F J 



(+) 



to be read from below upwards . The process must be initiated at p y N ' ^ 0. The 
vanishing Cl-\ = does not play any role in it. 

Roots of the SES secular determinants 

1.3. Dimension-independence. In every partial wave, the size L of the ma- 
trix in eq. 6JI coincides with the degree of the resulting secular polynomial. This 
value does not change with the growth of the dimension N + 1 of the vector p. 
Thus, the degree N of our polynomial wave functions enters our tridiagonal secular 
equation as a mere parameter. Containing extremely elementary matrix elements, 
the fully explicit form of this equation reads 

/ (L-l)b-F 1-L \ 



det 



-2N (L-3)b-F 
-2N + 2 



V 



4-2L 

(3 — L)b — F 
-2{N + 2- L) 



1 - L 
(1-L)b-F ) 



and specifies the family of the admissible charges F — F^\ at any integer L > 
and index k = 1, 2, . . . , L. The first few partial waves are exceptional since their 
eigencharges may in principle be defined by closed formulae a,t L — 1, L = 2, L = 3 
and L = 4. At the simplest choice of L = 1 (which may mean both the s— wave 
in three dimensions and an even state at D = 0), one does not obtain anything 
new. Secular equation provides the single root Fj^) which is equal to zero 

at all N. In the limit e — » 0, the quasi-even L — 1 solutions converge to the well 
known Hcrmitc polynomials. One just re-discovers the exactly solvable harmonic 
oscillator basis at even parity. We must re-emphasize that at any L > the number 
of our quasi-even roots does not change with the growth of the dimension N. This 
is of paramount importance since the practical determination of eigencharges is 
performed for all the indices N at once. Each of these families numbered by k 
contains infinitely many elements numbered by the first subscript. This might 
facilitate their future applications (cf. section 2 below). 

1.4. The first nontrivial generalization of oscillator basis: L = 2. Let 

us move to the index L = 2 giving the p— wave in two dimensions or to the s— wave in 
four dimensions. The related t — 1/2 wave functions may be chosen as compatible 
with the even-quasi-parity criterion H1.5|) . 

1 



p ( Q + \N,k) 



2N 



[b + F N , k ] p{ +) (N,k), 



1,2 



^we abbreviated /3„ 



- (2n + 1 — L) b for the sake of brevity 
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Secular equation [jl.fij) reads 

det ( -2^ _ b l F )=F 2 -b 2 -2N = 

and has two roots 

(1.7) f£> = ^W+2N), F$ = -^ 2 + 2A). 

[Note the contrast with the necessity of searching for roots of the polynomial (|4.4|) 
of the N—L+ 1st degree!] We encounter the non- vanishing charges for the first time. 
They grow with the increasing size of the shift b and with the quantum number N 
(i.e., with the energy). One finds their set more similar to the "F — line" of the 
harmonic oscillator than to the QES quasi-odd roots of Appendix B. Indeed, in the 
latter case the roots F^ ^ must be determined by the method which depends on N. 
In this sense we may now speak about the most natural and unique non-Hermitean 
generalization of the harmonic oscillator even-parity basis of £2(0, 00). 

1.5. Cardano charges at L = 3. At L — 3 the comparatively compact form 
of our secular equation 

/ 2b -F -2 \ 

(1.8) det -2N -F -2 =0 

\ -27V + 2 -2b -F J 

enables us to search for the triplet of charges {-F/v,i, Fn,2, Fn,3} via the closed (so 
called Cardano) formulae. One of alternative strategies may consist in a transition 
from the SES eigencharges F — F N _ k (b) to the inverse functions b = bjy,fc'(F). 
Such a trick lowers (by one) the degree of the secular polynomial at the odd values 
of L and, hence, extends the solvable class up to L = 5. A serious shortcoming 
of such an approach may be seen in the parallel deformation of the energies (|1.4|) 
which would change with the shift b. 

In the similar spirit, another simplification of formulae may be based on the 
formal elimination of N . The consequences may be illustrated by the secular L = 3 
determinant i|1.8fl which represents the mere linear problem for N, with the unique 
solution 




This formula suggests a simultaneous tuning of both the shifts ^[1.2] an< i the related 
charges F[ 12 ]- This may be achieved, say, by their hyperbolic re-parametrization 
with F(t) = V2e l and coshi = \ZjV cosha(t) etc. In terms of the parameters 
t = tN,k and/or a = ajv,fc, the spectrum of energies becomes deformed by the 
induced parametric dependence of b = bnwft) — ±v2 cosh 2 t-2A - e -t /V2. 
Purely formally, the eliminaton of N might extend the use of closed formulae up 
to L = 9. 

Higher partial waves 

1.6. Numerical methods. A shortcoming of the present SES construction 
lies in the growth of its complexity at the large L. At L > 5 one already cannot 
generate closed formulae for eigencharges, and the purely numerical search for the 
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roots F is necessary at all the very large L ^> 1. In an illustration using L — 3, and 
6 = 5 one gets the three eigencharges 

F = {10.757, -10.400, -0.35755} 

at the smallest possible TV = 2. They smoothly grow to the values 

(1.10) F = {89.98, -89.975, -0.0049407} 

evaluated at very large TV = 1000. This type of calculation is very quick and its 
results exhibit a smooth TV— dependence sampled in Table 2. 

Table 2. TV— dependence of eigencharges at L = 4 and 6 = 5 



TV 


Fn,u 


3 


-15.611 


-5.9279 


4. 888 7 


16.651 


30 


-27.149 


-9.290 9 


8.929 4 


27.511 


300 


-74.856 


-24.984 


24. 936 


74. 904 


3000 


-232.82 


-77.610 


77. 605 


232. 83 


30000 


-734. 99 


-245.00 


245. 00 


734. 99 



1.7. Large TV expansions. Any information concerning the eigenvalues F = 
F(b) and/or b = b(F) may shorten the necessary computations. At the intermediate 
L = 5 (which corresponds to the d— wave in three dimensions) we may eliminate, 
for example, the product 512 FN^ which is equal to 



-7686 - 256F6 2 + 768F + 40.F 3 ± 24^(10246 2 + 1926F 3 + 512F 2 + F 6 ) 

and use this formula as a constraint. Fortunately, closed formulae of this type may 
also open the door to the perturbative methods. 

Empirically, many numerically computed roots F^^ are very smooth functions 
of TV, especially in the asymptotic domain where they may be well approximated 
by their available large— TV estimates. For example, the exact result (|1.1U|1 already 
lies very close to its leading-order analytic estimate 

F « {V8TV, -V8TV, -6/TV} w {89.44, -89.44, -0.005}. 

At the intermediate values of the dimension TV, precision may be still insufficient 
but one can easily evaluate the large— TV corrections. For illustration, let us return 
to eq. i|1.9fl an d re-write it, in the case of its smallest root F = —bF/N , as the 
strictly equivalent formula 

^= 1+ TV^ 3 ' K + (6 2 6 1)/2TV ' L = 3 - 

It is suitable for iterations which represent our root as the following power series, 

8 3 2 3 d i 

TV 3 TV 6 TV 9 TV 12 
Also the other two L — 3 roots may be represented by the similar asymptotic series. 
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At the higher partial waves the same method works with a comparable effi- 
ciency. At L = 4, for example, the use of the variable M = y/2N + b 2 — 2 > 
compactifies secular equation <|1.6|l . 

(F 2 - M 2 ) (F 2 - 9A1 2 ) = 36- 48bF . 

In accord with Table 2, its asymptotically dominant 0(M 2 ) part determines the 
four distinct leading-order asymptotics of F ~ g M where g = ±1,±3. Once we 
re-normalize F = g M y/1 + R, the four exact charges obey the relation 

mVTTR 1 36 1_ 

~ 8|(?|-16 + (? 2 i? qM 3 + 8\g\- 16 + g 2 R g 2 M 4 ' ~ ' 
Its iterations evaluate the correction term R with an astounding efficiency. 



2. Quasi- VARIATIONAL SOLVABILITY on model spaces 
"Sufficiently large" finite subspaces in Hilberet space 

2.1. Left and right SES solutions. A full-fledged applicability of the stan- 
dard sets of QES wave functions is weakened by their incompleteness. Moreover, 
all the Hermitian QES Hamiltonians are usually interpreted as possessing just a 
finite multiplet of elementary bound states (at a given, special QES coupling) or 
Sturmians (cf. their example in Appendix B). For this reason, even the study of 
their own small perturbations is not easy at all jZ5| . Still, the sets of QES solutions 
themselves become large enough when we take into consideration both their fixed- 
coupling and fixed-energy subsets. In practice, such a philosophy proved useful for 
approximation purposes |BC| . 

In this section, we are going to return to the latter idea in the present broader, 
PT symmetric SES context. For the sake of dcfinitcncss let us re-consider Hamil- 
tonian of cq. as if it were defined at any charge F standing at the Coulomb-like 
force W(r) — ir~~ x . Then we may denote H — H{F) — H(0) + FW and, in the 
spirit of the Dirac's bra and ket notation, abbreviate |"0m i) = l-^i This enables 
us to re-write the special SES version of eq. in shorthand, 

(2.1) H(0) \M, k) + W \M, k) F M , k = \M, k) E M 

k = l,...,L, M = 0,1,... . 
Our PT symmetric Hamiltonian is non-Hermitean so that its left and right eigen- 
states will differ in general. One has to complement eq. (|2.1() by its counterpart 

(2.2) ({N,j\ H(0)+F NJ ((NJ\W = E N ((N,j\, 

j = l,...,L, N = 0,1,... 

where the same operators act to the left. After a return to the differential form of 
this problem i|l.lfl . a transition between Hermitian-conjugate equations (|2.1() and 
(|2.2|l may be re-interpreted as a certain reflection R in the space of the parameters, 

(2.3) R : {F < — ► —F, b < — ► -b, e < — ► -e .} . 

In an illustration using L — 2, our SES states of section 1 may be represented by 
a discrete set of points in the energy-charge plane. The values of their coordinates 
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(En,Fn^) are specified by the respective closed formulae ljl.4|l and l|1.7(l . These 
points may be perceived as located on the two (viz. left and right) branches of a 
single hyperbolic curve. These two infinitely large multiplets are mutually related 
by the transformation R of eq. l|2.3[) so that R\N, 1) = \N, 2)) and vice versa. 

2.2. Generalized Sturmian SES states as a basis. There is no a priori 
reason for a bi-orthogonality between our multi- indexed bras {{N, k\ = ((A\ and 
ket vectors \N',k') = \b) so that their overlaps 

Q A ,b = ((A\b) 

form a non-diagonal and asymmetric matrix in general. We have to assume that 
after a finite-dimensional truncation, this overlap matrix is invertible. Only in such 
a case we may define the inverse R = Q^ 1 and introduce an identity projector in a 
"sufficiently large" subspace of Hilbert space, 

(2.4) /= Yl \a)R a , B ((B\. 

aeJ kct ,BeJ bra 

In a constructive mood, let us return to equations Ij2.1|l and (|2.2(l and imagine that 
they share all their eigen-energies and eigen-charges. This enables us to write down 
the following two alternative matrix equations 

(2.5) ( (N, j | H(0) \M, k) = ((N,j\M, k) E M - ((N,j\ W \M, k) F M>k 

(2.6) ((N,j\ H(0) |M, k) = E N ((N,j\M, k) - F N>j ((N,j\ W \M, k) 
with {N,j) G Jbra and (M,k) € Jket- Their subtraction gives the constraint 

(2.7) (F Mtk - F Nd ) {(N,j\ W \M, k) = {E M ~ E N ) Q(iv,i),(M,fc) ■ 

Due to the way of its derivation, this relation may be understood as an immediate 
generalization of bi-orthogonality. In particular, we see that within the subspace of 
a single Sturmian multiplet (i.e., for M ~ N), the left-hand side expression must 
be a diagonal matrix with respect to its second indices. For our present purposes 
we shall abbreviate ((N,j\W\N,j) = wjvj. 

Non-QES bound-state problems 

2.3. Matrix equations. The knowledge of a basis can facilitate the study of 
many perturbed Hamiltonians via textbook perturbation recipes [F]. The use of 
the present SES bases might lead to a new progress in this area. Let us recall the 
respective right and left forms of any Schrodinger non-QES bound state problem 
with the above-mentioned structure, 

[H(0)+FW] |*> =E(F) |*) , 

((#| [H(0) + F W] = E(F) ((*| . 
Assuming that F ^ F^ ES ^ and using eq. 12.4fl we may insert 

|*) = ]T \a)R a , B {{B\V) = J2 \ a ) h « > 
oe JfcetjBe j bra a,eJk B t 



«*|= Yl ma)R a , B ((B\= 9b((B\ 

a£j k<lt ,BeJ bTa Be Jbra 
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and arrive at the double eigen-problem 

]T {(A\H{F)\b}R btC ((C\*) = E((A\V), A e J bra 

(plus its - omitted - conjugated companion), i.e., at the two conjugate linear alge- 
braic systems of equations written in terms of the same matrix, 

(2.8) Z(E,F)h = 0, g t Z(E,F)=Q 

Z Atb (E, F) = ((A\ H(0) \b) - E((A\b) + F ((A\ W \b) . 

Now, one could introduce a small parameter A = F — F^ ES%> and try to construct, 
say, E — E{\) in the form of a power series in A. 

2.4. Matrix elements. Long before any numerical determination of bound 
states, we must evaluate all the necessary matrix elements as an input. In practice, 
the latter step is usually the most time-consuming part of the algorithm. One has 
to optimize it in the present QES setting, therefore. 

In the first step we recall eq. <|2.6|) and eliminate all the matrix elements of 
H(0). This means that in eq. (|2.8|l we reduce the costly input to the mere evaluation 
of the matrix elements of the weakly singular Coulomb potential W(r) = i/r, 

Z Ajb (E, F) = (F- F A ) ((A\ W \b) - (E - E A ) ((A\b) . 

In the second step we keep M ^ N (i.e., we stay out of the Sturmian subspaces or 
diagonal blocks in the matrix Z) and postulate the absence of a random degeneracy 
of charges. This means Fm.u Fn,j so that we are permitted to re-arrange the 
bi-orthogonality-like relation 12.7(1 into a further reduction of the necessary input 
information, 

((TV, j | W |M, k) = J? M - En Q {N , ]U M,k) , M^N . 

In this way we arrive at the final form of our linear Schrodinger non-QES algebraic 
problem for the right eigenvectors, 

WN,j hN,j + ^2 ~~rp~~~~ Q(N,j),(Kffi) h K,p = 

K{^N),p J ' p 



_ E - E N ^ 

- p _p 1^ y(N,j),(M,k) h M,k , 

N >3 M,k 

j = l,2,...,L, N = 0,1,..- ■ 
For left eigenvectors, the system of equations is very similar though not equivalent, 

\ - E N - E K 

W N ,j 9N,j + 2^ ~P P Q(K,p),(N,j) gK,p = 

K(=iN),p J ' p 

E — En \ - „ 

- p _p 2^ y(M,k),(N,j) 9M,k , 

Nj M,k 

j = l,2,...,L, JV = 0,1,... . 

To solve any of these systems, say, by a perturbation technique, we just need to 
know the overlaps Q and the vector of Coulombic elements wn,j- 
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3. Concluding questions and remarks 

3.1. Does the quantum quasi-exact solvability have a classical ana- 
logue? In some studies a new relationship has been traced between classical and 
quantum mechanics. Its core may be seen in a correlation between the concepts 
of integrability and solvability. Needless to repeat, the former feature plays a key 
role in classical systems while its versions known as superintegrability and maximal 
superintegrability acquire more relevance in quantized world. 

For many quantum Hamiltonians the latter qualities also seem related to a 
purely algebraic property of quasi-exact solvability R . In such a context, it is ex- 
tremely exciting to ask questions about the robustness of the latter correlation and 
about the natural ways of the definition of the solvability itself. In a purely prag- 
matic manner, PT symmetrization could offer here another bridge towards the new 
solvable models. It is evident that after one weakens the traditional requirement of 
Hermiticity, many new QES models may be found, indeed [CIl IZ6| . 

In the present paper we have seen that for the D dimensional central Coulomb 
plus shifted harmonic oscillator its quasi-exact solvability acquires a fairly unusual 
modified form. A significant difference appears between the quasi-even and quasi- 
odd states, where only the properties of the latter set remain standard. The facili- 
tated construction and unexpectedly non-numerical character of the former quasi- 
even family make it similar to the current complete oscillator basis on half-axis. 
This locates our new "semi-exact" solutions in a gap between their older exact and 
quasi-exact neighbors. 

3.2. Do we need more semi-exactly solvable models? Our presentation 
of details started from the reduction of the four-term recurrences of our "paper 
I" |Z4j to their improved three-term form. This was rendered possible by our new 
ansatz which also proved able to reproduce all the available older results. We under- 
lined that one of the most characteristic shortcomings of the standard QES equa- 
tions lies in the necessity of solving the adjacent linear algebraic N— dimensional 
eigenvalue problem of growing size N which, in our particular example, selects the 
admissible charges. 

In this context it is important that we succeeded in revealing the existence 
of a new, "facilitated" or "semi-exact" elementary solvability emerging when the 
quasi parity was assumed even. Infinitely many of these states acquire the exact 
polynomial form in each (= to— th) partial wave and at any element E = E(N) of 
an equidistant set of the energies. In fact, our states form the Sturmian L— plets 
(numbered by N = 0, 1, . . .) at charges i*jv,fc numbered by k = 1, 2, ... , L where 
L = D — 2 + 2to is N— independent. 

The well known even-parity basis on £2(0,00) (made of Hcrmitc polynomials) 
is re-obtained in the simplest L — 1 special case where the resulting (single) QES 
charge is zero, F = 0. In the first nontrivial L = 2 case we have got the two 
alternative eligible infinite series of states at charges F — zty/b 2 + 2N (where b is 
the shift). We have shown that and how the similar sets might serve as a source of 
the new matrix reformulations and approximation techniques within perturbation 
theory or variational considerations. 

Our main attention has been paid to the generalized QES ansatz admitting 
solutions with both quasi-parities. Our main result, viz., a new version of the SES 



as reviewed, e.g., by Pavel Winternitz in this volume IWil 



12 



MILOSLAV ZNOJIL 



recurrent construction is "almost exact", first of all, due to its N— independence. 
In the first few lowest partial waves our infinite families of semi-exact states proved 
very transparent and close, by their applicability, to the basis of harmonic oscilla- 
tors. We revealed that their construction is amenable to an efficient approximative 
treatment even at the higher partial waves via the Taylor-series expansions simi- 
lar to the common large— N perturbation theory. This might not only initiate the 
study of the old models in the non-QES domains of couplings but also offer a new 
motivation for an intensified search for the new SES models. 

3.3. Will our new SES bases find efficient practical applications? The 

matrices of overlaps QA,b are, presumably, non-diagonal. In the other words, our 
SES basis states \N, k) and \N,k)) are not mutually bi-orthogonal. Any deeper 
insight in their overlaps would be appreciated in applications, therefore. In the 
future, this could facilitate, say, a search for a non-Hermitian analogue of the QES 
+ extrapolation trick as suggested by Burrows et al |BC[ for the Hermitian sextic 
anharmonicities . 

Our formalism is not yet fully prepared for initiation of any practical numerical 
or perturbative calculations in the non-QES regime of course. At the same time, 
we already succeeded in a drastic reduction of the number of the necessary input 
matrix elements. This supports our belief that our present constructions represent 
a valuable key step towards deeper understanding of relations between the degree 
of integrability in classical mechanics and an extent of exact solvability in quantum 
mechanics. At this moment we have to re-emphasize the definite progress in our un- 
derstanding of properties of the "forgotten" quasi-even solutions, their quasi-exact 
solvability of which may be characterized as "significantly facilitated" . They are 
distributed in the coupling-energy plane along curves which may be interpreted as 
a common generalization of the usual bound-state straight lines 5 and the straight 
lines of the so called Sturmian states 6 . In this sense, their families may be under- 
stood as certain generalized Sturmian families defined on some curves in the plane 
of charge and energy. These sets contain infinitely many elements and definitely 
resemble and generalize the current types of bases. 

3.4. Could the complexification serve as a regularization recipe for 
SW models? In connection with the strongly singular behaviour of many classi- 
cally superintegrable models at r = 0, a new role of the present weakening of their 
Hermiticity may be also sought in its parallel regularization effect. In this role, 
the PT symmetric regularization proved very useful and really vital in the vari- 
ous PT symmetrized versions of the so called supersymmetric (SUSY) quantum 
mechanics |ZCl IE] . For example, in accord with refs. |Z8| the SW Hamiltonians 
H — p 2 + r 2 + G/r 2 may form the formal SUSY partners H^ L ^ R y Once their sin- 
gularity at r — is circumvented via the complex shift of the coordinate axis, their 
spectra remain real and discrete and exhibit still the usual SUSY-isospectrality 
pattern in the so called unbroken PT symmetry domain with G > —1/4. 

In Z9] we have shown that beyond the latter domain, the PT symmetry it- 
self is completely broken but the spectra stay partially real on a SUSY-generated 



^characterized by the fixed couplings and variable energies and exemplified by the spectrum 
of the exactly solvable harmonic oscillator or by the Singh's sextic QES oscillators |SB1 

^characterized by the variable couplings and fixed energies and exemplified by the exactly 
solvable zero-order models in ref. IZ7I or by the Flessas' non-polynomial QES oscillators IWFI 
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set of complex G. In addition, an alternative SUSY scheme with equal spikes, 
G(l) = G(R) may be introduced giving an innovated super-Hamiltonian factorized 
in terms of certain creation and annihilation differential operators of the second 
order which conserve the quasi-parity and are mutually not adjoint. Together with 
the original Hamiltonian H the latter two operators commute in accord with the 
Lie algebra sl(2, R). 

All these preliminary results and hints are encouraging and suggest that in the 
context of the Calogero and/or Smorodinsky-Winternitz supcrintegrable models 
with strong barriers the PT symmetry-induced tunneling effects did not say their 
last word yet. 
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4. Appendices 

A. Charged harmonic oscillator and the concept of quasi-parity 

As we already mentioned, the shifted and charged harmonic oscillator of ref. 
[H] is a characteristic illustrative example of QES model in quantum mechanics. 
Its Hcrmitian Schrodinger equation 



(4.1) 



d 2 | 1(1+1) | / 
dr 2 r 2 r 



+ 2 gr + r 



ip(r) = E (p(r), r € (0, oo) 



has been replaced by its FT— symmetrized version (ll.lll in ref. |Z4| . After the 
severe restriction of the potential to s— waves in three dimensions (i.e., to D = 3 
and I = 0, motivated by technical reasons) we succeeded in confirming its quasi- 
exact solvability in non-Hcrmitian regime. In the present continuation of paper 
|Z4| this result is extended to all the partial waves I = (D — 3)/2, (D — 3)/2+ 1, . . . 

and dimensions D = 3, 4, 

The first hint indicating the open possibilities in connection with eq. (|4.1() 
appeared in ref. |Z2j where the concept of quasi-parity has been introduced. One 
was able to conclude there that all the solutions of the differential Schrodinger 
equations of the radial type prove separated into quasi-even or quasi-odd states 
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( x ) m a wa y which complements the quasi-odd solutions of ref. |Z4| by their 
present "forgotten" partners with quasi-even symmetry. The key notion of the 
quasi-parity itself may be defined either via asymptotics or, more easily, via the 
specific behaviour of the wave functions near the singularity at r[x) ~ 0, 

ip£~'(x) ~ (x — ie)~ l , i>^\x) ~ (x — ie) e+1 , x^ie. 

The name suggests that the quasi-parity coincides with the ordinary parity after 
limiting transition e — » to a Hermitian Hamiltonian. In this limit, the quasi-even 
solutions themselves need not lose their normalizability (and may remain physical) 
whenever the strength of their singular repulsion is sufficiently small (cf. an explicit 
illustration in ref. |Z8p . 



B. QES states with odd quasi-parity 

Formula l|1.4|) for the energies is independent of the charge and was already 
known to Hautot H . Thus, the energy is a constant for all the solutions of eq. 
(|1.3() . The multiplets of bound states of such a type 7 find applications, say, in 
perturbation theory |Z7j . Up to a free normalization pn ^ 0, their present closed 
definition 
(4.2) 



/ B 
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B N -1 

A N 



Cjv-i 
B N 



is unique at any energy E, shift b and charge F. 

Our new ansatz i|1.2fl is sufficiently flexible and reproduces all the older quasi- 
odd QES solutions of refs. [HI and |Z4| . Indeed, we may mark them by the 
superscript and characterize their set by the following boundary condition, 



(4.3) p%->=p[->=... = = 

which reduces the range of the indices in formula l|4.2|l . j = N, N — 1, . . . , L . We 
have to guarantee that the QES recurrent recipe terminates, i.e., that the value of 
p^ L \ vanishes. This condition has the secular form 



(4.4) 



= det 
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Cl+i 
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N-l 



B N -x 
A N 



Cn-i 
Bn 



and its occurrence is characteristic for all QES systems [U]. At every admissible 
energy it specifies a multiplet of the admissible charges F = F^ 'j numbered by the 
second subscript j = 1,2, . . . ,N — L + l. At each N and L our secular polynomial 
is of degree N — L+l. Of course, the practical determination of the quasi-odd QES 



^defined at TV different charges i*jv,j and called Sturmians 
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eigencharges is a purely numerical procedure for N > L + 3. Fortunately, as we 
have seen in section 1, this shortcoming may be suppressed in the quasi-even case. 
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